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Abstract 

We compute the multipartite entanglement measures such as the global entanglement of various 
one- and two-dimensional quantum systems to probe the quantum criticality based on the matrix 
and tensor product states (MPSs/TPSs). We use infinite time-evolving block decimation (iTEBD) 
method to find the ground states numerically in the form of MPSs/TPSs, and then evaluate 
their entanglement measures by the method of tensor renormalization group (TRG). We find 
these entanglement measures can characterize the quantum phase transitions by their derivative 
discontinuity right at the critical points in all models considered here. We also comment on the 
scaling behaviors of the entanglement measures by the ideas of quantum state renormalization 
group transformations. 
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I. INTRODUCTION 



Quantum phase transitions (QPTs) are conjectured to happen in many interesting physi- 
cal systems and attract much attention in modern condensed matter physics |l| . The QPTs 
occur at zero temperature when a parameter is tuned to a critical value, and are facilitated 
not by the thermal but the quantum fluctuations in the ground state. Especially, the correla- 
tion length diverges at the critical point and some order parameters characterize the change 
of the ground states across the critical point. The diverging correlation length also implies 
the different parts of the system are strongly correlated quantum mechanically. Therefore, 
the quantum correlation and entanglement play a key role near quantum critical point. In 
other word, the ground state is highly entangled and the QPT should be characterized by 
the abrupt change of the quantum entanglement 2j, [3] , mainly probed by the single-site von 
Neuman entropy (1-tangle) or the two site concurrence for one-dimensional (ID) quantum 
spin systems, for a review see j^J]. 

On the other hand, there are very few multipartite entanglement measures which can be 
used to characterize the QPTs, among which the global entanglement (GE) was proposed to 
probe the quantum criticality and has been tested for ID quantum spin systems {5, 6, 7, §}. 
However, as far as we know, almost nothing has been done to use GE to probe the quantum 
criticality of the 2D quantum spin systems, except a related fidelity computation for trans- 
verse Ising model in 9fl. This is mainly due to the lack of the exactly solvable or numerically 
tractable ground states in 2D case. The situation, however, changes recently because of 
the new development of the generalization of the ID density-matrix renormalization group 
(DMRG) 10] and matrix product states (MPSs) U, 12, Q, Q to the 2D tensor product 
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states (TPSs) [15|, 

In this paper, we would like to probe the ID and 2D quantum criticality numerically by the 
multipartite entanglement measure such as GE based on the approximation of MPSs/TPSs 
to the true ground states of the systems. Mainly, we will utilize the infinite time evolving 
block decimation (iTEBD) [ijl numerical method to find the approximate ground states 
in the form of MPSs/TPSs, and then use tensor-network renormalization group (TRG) 
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method to evaluate the entanglement measures of the corresponding TPSs. By exploiting 
the translational invariance of the ground states, the TPS and TRG algorithm become quite 
numerically tractable 17. ["181 1 1 



Below we will briefly summarize the entanglement measures we will evaluate in this paper. 
The first is GE which is geometrically defined as the distance between lib) and the nearest 

nn 

separable state \cf>) [5|, |6|, i.e., 

Ety) = -\ogk 2 max , A max = max (1) 

where the maximization of the fidelity A is with respect to the variation of In the 
representation of TPSs, the fidelity can be computed by using TRG. Next, the TPSs with 
translational invariance can also be seen as the bipartite Schmidt decomposition of the 
ground states, whose Schmidt number Aj's (i = l,...,d s with d s the physical dimension of 
the spin.) are the singular values obtained from the singular value decomposition (SVD) 
of the site tensors in performing the iTEBD. Then, we can straightforwardly evaluate the 
bipartite entanglement measure 

^P = -^^logA, 2 . (2) 

i 

Finally, for comparison we will also evaluate the single-site von Neumann entropy (1-tangle), 
which is denoted as S±. 

Our results show that all three above entanglement measures are equally good to char- 
acterize the quantum criticality by their derivative discontinuities right at the critical point 
for various ID and 2D spin systems considered in this paper. This helps to establish the 
universal nature of entanglement probe for quantum criticality. Besides, we also study the 
scaling behaviors of the entanglement measures for ID and 2D systems based on the idea of 
quantum state renormalization group (RG) transformation. In ID case we find the expected 



scaling behaviors as done for different models in [36|, |37j, |38|, however, the method fails for 
2D case. 

This paper is organized as follows. In the next section we will evaluate the above three 
entanglement measures for various ID spin systems based on MPSs formalism, including 
spin 1/2 XY model and spin 1 XXZ model. In section III we will evaluate the entanglement 
measures for various 2D spin systems by the TPSs and TRG algorithm, including transverse 
Ising model, XYX model and XXZ model. In section IV we discuss the scaling behaviors 



of the entanglement measures for both ID and 2D systems based on quantum state RG 
transformation method. We then briefly conclude in section V. 



II. ENTANGLEMENT MEASURE IN ID SPIN SYSTEM 

For ID quantum many-body systems, the ground states could be expressed in terms of a 
matrix product state (MPS) as 

|V>>= £ Tr[A^% 1 )A^(s 2 )...A^(s N )}\s h ...,s n ) (3) 

S1,S2...,SJV 

where Sj = l,...,d s for i = 1,...,N, and A^(si) are Xi-i by Xi matrices, with d s denoting 
the physical dimension and Xi denoting the dimension of the z-th bond. 

We use the iTEBD method to numerically solve for the above MPS, namely, by acting 
the imaginary time evolution operator exp(— tH) on an initial state \ipo). The operator 
exp(— tH) for small enough r can be expanded through a Suzuki- Trotter decomposition 
as a sequence of two-site gates U^ l ' l+1 ^. Assume translational invariance the MPS can be 
formally written in the following form \%p) = Yl Sl s 2 s N TV[r^ i A j4 r^A' B ■ • s n ), where 

Tj's are x by x matrices, and A's are x by x diagonal matrices of singular values. Then, we 
only need to update four tensors T A , T B , X A , and X B for each evolution step. In the limit 
r — > oo, exp(— TH)\ip Q ) will converge to the ground state of H we are solving for. We also 
compute the magnetization M x = (ipg\crf \ip g ) and the global entanglement per site by the 
same method. 

The global entanglement can be computed according to ([T]), first we should maximize the 
fidelity (CQ) between the quantum state \ip) and a separable state For a system of spin 
1/2, the separable state |0) = <S)f =1 (cos(9i)\0) + sin(6 , i )|l)) where |0) and |1) are eigenstates 
of a z . So, the fidelity takes the form 

|(V#>| = \Tr(Tg^W 2] -Tg [N] )\ (4) 

where Tg® = E Si ^ [i] (^) ® B^fa) are transfer matrices, with B®(0) = cos(6>i), B®(1) = 
sin(6 l j), and A^(si) given in Eqns.Q. Again, assume translational invariance, 9i := 9 
and := Tg, then the fidelity can be reduced to A(0) = \Tr((Tg) N )\ = |S (A a ) iV |, 
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where A a are eigenvalues of transfer matrix Tg. In case the ground state is non-normalized, 

we should normalize the fidelity as li^Mj where 1(^1^)1 = \Tr(Tn^Tn^ ...Tn^)\ with 

yi(V#)l 

the transfer matrices Tnfi = ~52 s . A^(si) (g) A^(si), or assuming translational invariance 
|(^|^)| = \Tr(Tn N )\ = \H e (X e ) N \ with A e 's the eigenvalues of Tn. 

Finally, the maximization of the fidelity A(9) with respect to the angle 9 can be computed 
by numerical algorithms, and the global entanglement follows, i.e., 



IV (\ \N\2 

£W = - log S = - log TO- (5) 

A. Spin 1/2 chain 

We consider the quantum phase transition of a spin 1/2 quantum ID XY closed spin chain 
with transverse magnetic field. The Hamiltonian is (setting the ferromagnetic coupling to 
one.) 

n n 

H = - E 2 [(1 + l)a ^ + (1 - 7)(T ^+ i] + h J2°i ( 6 ) 

where a", a = x,y,z are the Pauli matrices acting on site i, and h > is the strength of the 
transverse field. The model reduces to the quantum Ising model for 7 = 1. We will focus 
only on the parameter range < 7 < 1. The ground state for the XY spin chain can be 
solved analytically, which exhibits three phases: the oscillatory (O), the ferromagnetic (F) 



and the paramagnetic (P) ones, and the phase changes go as follows 41[. As we turn on 



and increase h, at h = y 1 — 7 2 the ground state will change from O to F phase if 7 7^ 0. 
This, however cannot be captured by the order parameter M x because it is a cross-over not 
a phase transition. Increase h further up to h « 1, we will see the transition from F(or O 
if 7 = 0) phase with nonzero M x to P one with zero M x , which is a second-order quantum 
phase transition in the thermodynamic limit. 

Here we use iTEBD method to compute the ground state wave function in the form of 
MPS for system size N = 2 10 , and the result is plotted in Fig. 1. Besides, in Fig. 2, we plot 
the entanglement measure to characterize the quantum phase transitions, which reproduces 

n 

the analytical results in |6|. 

In Fig. 2(a) we note that near the critical point the GE has a maximum value which 
implies stronger correlation near the critical point, moreover, the derivative of GE with 
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FIG. 1: The magnetizations v.s. transverse magnetic field of XY spin chain using MPS with bond 
dimension \ = 16. It indicates a 2nd-order phase transition from F to P phase at h c ~ 1.005(Ising, 
7 = 0) and h c — 1.005(XY, 7 = 1/2). However, it does not capture the phase transition from O to 
F phase. 
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FIG. 2: Entanglement measures v.s. transverse magnetic field of XY spin chain using MPS with 
bond dimension x = 16- (a) (Left) GE v.s. h. Note that near the critical point h c ~ 1.005, the 
first derivative of GE shows log \ h — h c \ scaling behavior, see the boxed graph. This agrees with 

n 

the observation in [6]. (b) (Right) Sbp v.s. h. 

respect to h is discontinuous and diverges at h ~ 1.005. The discontinuity of a physical 
quantity such as free energy or its derivatives is the notion of phase transition. Here the 
derivative discontinuity of the entanglement measure plays the similar role in characterizing 
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the phase transition. On the other hand, for 7 = 1/2, at h 2 + 7 2 = 1 the entanglement 
density vanishes and its derivative is discontinuous, it then characterizes the oscillatory (O) 
to ferromagnetic (F) phase. In Fig. 2(b) the bipartite entanglement entropy Sbp shows the 
same behaviors. Moreover, Sbp is far easier to evaluate than GE since we only need the 
the singular values of MPS for the ground state without further numerical maneuver for 
calculating the expectation values. 



B. Spin 1 chain 

To demonstrate more on the power of GE in characterizing the quantum phase transition, 
we now consider a spin 1 XXZ Heisenberg chain with anisotropy. The Hamiltonian is 

n n 

H = + S " S li + - WS^] + D Sf (7) 

i=i i=i 

where Sf, a = x,y,z are spin-1 operators. The parameter D represents the uniaxial 
anisotropy. The ground state phase diagram of this model consists of six different phases. 
Here we focus our attention on phase transitions between large- D phase, Neel phase and 



Haldane phase (characterized by nonzero string order parameters) [2l|, |22j, |23| . We consider 
two cases of quantum phase transitions: (1) transition between large-/} and Neel phases by 
tuning D but fixing J z ; (2)transition between Neel and Haldane phases by tuning J z but 
fixing D. We can numerically solve the ground state of the spin chain in the form of MPS 
for different values of D or J z . 

Obviously, the ideal large- .D state is 1 000 . . .), and the ideal Neel state is |1, —1, 1, — 1 . . .) 
or| — 1,1,— 1,1...). By increasing the uniaxial anisotropy D the system undergoes a first- 
order quantum phase transition from Neel to large-/} phase characterized by the staggered 
magnetization M s = -p J2iLi(~ ^Y(Sf)- Using the solved MPS, we plot in Fig. 3(a) the M s 
versus D with fixed J z = 2.59 for bond dimension x — 4 with system size N = 2 10 , and in 
Fig. 3(b) the M s versus J z with fixed D = 0.3. 

On the other hand, we plot the various entanglement measures in Fig. 4, which are also 
capable to characterize the quantum phase transitions mentioned above. To compute the 
GE for a translational invariant spin 1 chain with Neel order, we should choose the following 
separable state ansatz |0) = (\4>i) ® |02)) lX> ^, where \<fii) = (sin^ cos5»|l) + sin^ sin5j|0) + 
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FIG. 3: Staggered magnetization of spin 1 XXZ chain, (a) (Left) M s v.s. D with fixed J z = 2.59, 
which indicates a phase transition from Neel to large- D phase at D ~ 2.3. (b) (Right) M s v.s. J z 
with fixed D = 0.3, which indicates a phase transition from Haldane to Neel phase at J z ~ 1.17. 



cos^l — 1)) with |1), |0) and | — 1) are the eigenstates of the S z . We can then compute 
the global entanglement as usual. For comparison we also calculate the 1-tangle Si and 
the bipartite entanglement Sbp- Again, the derivative discontinuities of the entanglement 
measures characterize the quantum phase transition. Moreover, we see that all three en- 
tanglement measures do equally good job for showing the phase transition as the staggered 
magnetization does. 



III. ENTANGLEMENT MEASURE IN 2D SPIN SYSTEMS 

We now turn to investigate the global entanglement measure in 2D quantum spin system. 
Especially, we would like to see if it is capable of characterizing the quantum phase transition 
or not. 

As in ID case, we will compute the ground state wave function in the form of the tensor 
product state (TPS). Several methods have been developed in the past few years, including 
the variational optimization of the expectation value of the energy [3] and the iTEBD 



method 181 ] . The ground states of 2D lattice systems obtained by using any of the above 
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FIG. 4: The entanglement measures of spin 1 XXZ chain. Their derivative discontinuity charac- 
terizes the quantum phase transition as good as the staggered magnetization does, (a) (Left) Neel 
to large-!? phase; (b) (Right) Haldane to Neel phase. 



methods, can be expressed in terms of a TPS as 



tTr[A^(s 1A )A [1 ' 2] (si,2) ■ ■ ■ A^ N \s N , N )]\ Sl>1 , s 1>2 , s N , N ) (8) 



Sl,l,Sl,2,-, s JV,iV 



where the values of the physical spin = 1, • ■ ■ d s for i,j = 1,...,N ; tTr is to sum 
over all indices of tensors; and A\ ^\ M (s;.j)'s are rank-five tensors with the bond indices 
r,d,l,u = 1, • • • , X- We again call x the bond dimension and d s the physical dimension. 
In the following, we will use the iTEBD method to calculate the TPS on square lattice, 
which is quite similar to the ID case though we should update the tensors and the Schmidt 
values by more subtle singular value decompositions, for some example see [2^|. In order to 
minimize the Trotter error of iTEBD, we usually take 8r = 0.1. 

To compute GE, we need the separable state ansatz as |0) = <S^=i( cos (^,j)|0) + 
sin(6ij)\l)). Then, the norm of a TPS and its overlap with a separable state are 
= \tTr(Tn^ . ..Tn^ N ^) \ and |(^#)| = \tTr(Tg^ . . . Tg^ N ' N ^) \ respectively, where 



Tn 



E A 



1,31 (V 



B i s i,j) are the transfer matrices, and B^(0) := cos(9 i; j), B^ ,: >'(1) 



sin(6 l jj). For a 



translationally invariant system the fidelity A 



A 



a/KV#)I 

\tTr(Tg ®Tg®Tg& 
^J\tTr(Tn ® Tn ® Tn 



can be expressed as 



(9) 
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where Tg^l = Tg, Tn™ = Tn, Wi,j. We can then obtain the GE by vary it with respect 
to 9ij. 

Unlike the ID case where the matrix trace in (j3J) is computationally tractable. In 2D, the 
tensor trace tTr over all connected indices of the above norm or overlap is exponentially hard. 
Instead, a method called tensor renormalization group (TRG) is developed in [l7, I8.fi9l.l20l 
to make the double tensor trace polynomially calculable by merging the sites and truncating 
the bond dimension of the merging lattices according to the relevance of the components 
in the Schmidt decomposition. The cutoff of the merging bond dimension is denoted as 
D cut which controls the accuracy of the computation. Each step of the TRG will reduce the 
number of sites by half. Eventually, the tensor-network is reduced to only four sites and 
the double tensor trace for the norm or overlap can be calculated easily. For more detailed 
description of the TRG method, please see I7.li9l.l20l]. We will then adopt the TRG method 
to numerically evaluate the entanglement measure such as GE for the ground state in the 
form of TPS on the 2D square spin lattice. 



A. The 2D transverse Ising model 

To calibrate the accuracy of our numerical codes, we first consider the 2D spin 1/2 
transverse Is ing model on a square lattice, which is well-studied using TPS and TRG, for 



example see 19(. As a by-product we find that the entanglement measure such as GE can 



also characterize the quantum phase transition of the model. The Hamiltonian is 

<i,j> i 

where < i,j > stands for the nearest neighboring sites. 

By tuning h, there is a quantum phase transition characterized by the ^-symmetry 
order parameter M z = (ip\af\ip), and its occurs at h ~ 3.04 obtained from the unbiased 
quantum Monte Carlo simulation [24]. Our result using TPS and TRG is plotted in Fig. 
7. We consider the bond dimension \ = 4 for a system size 2 5 x 2 5 , and keep D cut = 16 
to ensure the accuracy of the TRG calculation. We find a second order phase transition at 
h ~ 3.25, which can be characterized either by the order parameter or the entanglement 
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FIG. 5: The 2D transverse Ising model: M z , GE, Sep and S\ v.s. the transverse magnetic field h 
for a system of size 2 5 x 2 5 . We adopt TPS and TRG with \ = 4 and D cu t = 16. The derivatives 
of these quantities all have a discontinuity at h ~ 3.25. In this figure, the GE and Sep are scaled 
up by a factor of 25 and 10, respectively. 

measures. Note that we obtain the bipartite entanglement per length Sep by averaging 
— Yl ^i,j 2 l°g ^i,j 2 over eight bonds of a plaquette. 



B. The 2D XYX model 



We now consider the 2D anti- ferromagnetic spin 1/2 XYX model in an uniform z-axis 
external magnetic field h, i.e., 



<ij> 



For A y < 1 the system has easy-plane (EP 



x-axis magnetization M x by tuning h 



behavior. The order parameter in the EP is the 



26|. 



The model shows two quantum critical points, and only one can be captured by the order 
parameter M x . This is similar to the phase diagram of ID XY quantum spin chain. Before 
turning on h, the ground state is in Neel phase with nonzero M x . As h increases to some 



critical value hf called factorizing field value (its theoretical value is hf — 2^/2(1 + A 



y) 



26|), the ground state becomes separable product state but still with nonzero M x . Therefore, 



the order parameter cannot characterize the cross-over but the entanglement measures do. 
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Increasing h further to the another critical value h c , there is a second-order phase transition 
from Neel to the disorder phase with zero M x . We also find that the entanglement measure 
can characterize this quantum critical point by their derivative discontinuities. Finally, as 
h —>■ oo, the system is in a polarized state with the magnetization of z-direction saturated 
to the unity. 

As in ID XY spin chain case, there is again a subtlety in choosing the separable state 
ansatz for evaluating GE in the Neel phase. By Neel and translational symmetry, we should 
parameterize the separable states via four different states, i.e., \<f>) = (\<f>i) <E> \<f>2) <8> {fa) <8> 
|</>4))® T with = cos(0j)|O) + sin(6'i)|l). We then compute the ground state and its 
entanglement measures by using iTEBD for TPS and TRG with \ = 4 and D cut = 16. 

We present our numerical results in Fig. 6 with A y = 0.25 and systems size of 2 5 x 2 5 , 
which shows a second-order quantum critical point at h c ~ 3.5 characterized by both the 
order parameter and the entanglement measures, and also the factorizing field value at 
hf ~ 3.16 at which the entanglement measures vanish as expected. Our results agree with 



the previous one based on 1-tangle 



The 2D XXZ model 



26] 



Finally, we consider the 2D spin 1/2 XZX model in an uniform z-axis external magnetic 
field h, i.e., 

H = -Y, w + - A(T M) - h E °f ■ ( 12 ) 

<i,j> i 



In the large spin-anisotropy limit (A > 1), it has been shown that [27|, |28|, |29( as h being 
increased, a first-order spin-flop quantum phase transition from Neel to spin-flipping phase 
occurs at some critical field h c . As we further increase h to another critical value h s = 
2(1 + A), the fully polarized state is reached. 

Here we consider the case with A = 1.5. Our results are plotted in Fig. 7 from solving 
ground state via iTEBD and TPS, and evaluating the order parameter and the entanglement 
measures. In Fig. 7(a) we plot the expectation values of the z-direction staggered magne- 
tization Ml, the uniform one M z , and the x-direction uniform magnetization M x . It shows 
that the Neel order M s z , will suddenly drop to zero at h c ~ 1.8 and remain zero after that. 
On the other hand, M z is zero for h < h c but suddenly start to grow at h — h c , and finally 
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FIG. 6: The 2D XYX model with A y = 0.25: M x , M z , GE, S B p and Si v.s. the transverse 
magnetic field h for a system of size 2 5 x 2 5 . We adopt TPS and TRG with x = 4 and D cu t = 16. 
The derivative of these quantities shows a discontinuity at h ~ 3.5, and the entanglement measures 
vanish at the factorizing field value h ~ 3.16. In this figure, the GE and Sbp are scaled up by a 
factor of 25 and 10, respectively. 

reaches a saturated value at h s ~ 5. Similarly, M x is nonzero only for h c < h < h s . Both M x 
and M z can be used to characterize the quantum phase transition from spin-flipping to the 
fully polarized state. Therefore, not all three order parameters can be used to characterize 
both quantum critical points. On the other hand, in Fig. 7(b) we can use either 1-tangle, 
GE or Sbp to characterize both quantum critical points by their derivative discontinuities. 
This suggests that all the entanglement measures discussed above may belong to the same 
universal class. This seemingly universal feature can be seen as an advantage to use en- 
tanglement measure as an index for quantum phase transition. Finally, unlike the usual 
derivative discontinuity of the entanglement measure for the phase transitions encountered 
before, we note that the first order spin-flop transition is characterized by the more dramatic 
change, namely, by the sudden jump of both order parameter and entanglement measures. 
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FIG. 7: The 2D XXZ model: (a) (Left) Magnetizations v.s. the transverse magnetic field h. The 
derivative of all three magnetizations has a discontinuity at h c ~ 1.8 for spin-flop phase transition, 
and only the M x and M z show discontinuity at h s ~ 5. (b) (Right) Entanglement measures v.s. 
h. Their derivative discontinuities characterize both quantum critical points. Here we adopt TPS 
and TRG with \ = 4 and D cu t = 16. 



This indicates that the degree of discontinuity of the entanglement measures as well as order 
parameter may be related to the order of phase transition. 



IV. COMMENT ON SCALING BEHAVIOR OF ENTANGLEMENT MEASURE 



In the usual renormalization group (RG) sense as the coarse graining, the quantum critical 
points are usually the IR fixed points of the original spin system, and around which the 
correlation length diverges and some scaling symmetries emerge. Theses scaling symmetries 
then dictate the scaling law of some physical quantities. On the other hand, it was known 
that the von Neumann entanglement entropy obeys the area law 35(, and especially in ID 
the entanglement entropy obeys the log scaling of the block size at critical point, with the 



coefficient proportional to the central charge of the under 



(CFT), otherwise it will saturate to a constant value 



36 



ying 1 + 1 conformal field theory 



37) 
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We may then wonder if the more quantum information motivated entanglement measures 
such as GE will obey the similar scaling law or area law or not. However, the straightfor- 
ward way of computing block entanglement measure is quite difficult because the physical 
dimensions of the block will grow with the block size exponentially, and make the numerical 
calculation formidable. Fortunately, in the framework of MPS or TPS one can adopt the 
method of quantum state RG transformation to merge the lattice sites [30[], by which the 
block entanglement measure of size L is the same as the single-site entanglement measure 



after L steps of quantum state RG transformation 38|, |39( . Moreover, at each step of quan- 
tum state RG transformation, the physical dimension is bounded by x 2 f° r ID and x f° r 
2D. However, this kind of truncation of physical dimension only throws away the local in- 
formations but still keeps the non-local ones intact. In this way, we can compute the block 
entanglement measure for MPSs/TPSs. 

In the following we will comment on the scaling behaviors of entanglement measure of 
MPSs/TPSs in ID and 2D quantum spin systems. 

A. ID scaling via MPSs 

When discussing the scaling behaviors of the block entanglement measure near the quan- 
tum critical point based on MPSs/TPSs approximation to the true ground state of the spin 
systems, we should be aware that the bond dimension % will play the role of the effective 
correlation length £ x of the MPSs/TPSs, which will not diverge at the critical point as for 
the correlation length of the true ground state £. The difference between £ x and £ is because 
we approximate the true ground state by the MPSs/TPSs with finite x- Therefore, £ x should 
be related to the bond dimension \. Indeed, in 
law for the ID Ising model, 



the authors found the following scaling 



tx ^ X 2 , (13) 
and the block von Neuman entropy (1-tangle) will be bounded by 

c c 

S x - q lo g£x - 3 lo SX (14) 

herefore, the typical scaling law 



351 ] will saturate to the value of 



where c is the central charge of the corresponding IR CFT. 
of the block von Neumann entanglement Sjj ~ ~ log L 31 
( THj) as L > £ x ~ x 2 ■ On the other hand, when away from the critical point, the correlation 

15 



length £ could be far smaller than £ x and the MPSs/TPSs approximate the true ground 
state well. 
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FIG. 8: The block entanglement entropy Sl for spin 1 XXZ chain with the fixed J z = 2.59 and the 
quantum critical point at D ~ 2.3. We adopt the MPSs of x = 2, 4, 8 for the ground state. (Left) 
Off the critical point, the block entanglement entropy saturates around L = 4. (Right) Very near 
the critical point, the block entanglement entropy saturates around L ~ x 2 , which can be seen 
more clearly from the numerical data. 



Based on the above discussion in the framework of the MPS/TPS, one may wonder if 
the scaling law (fl~3l) and the entropy bound (fT4"l) are universal or not for different ID spin 
systems and various entanglement measures. In Fig. 8 we plot the block von Neumann 
entanglement entropy of size L for the ID spin 1 XXZ chain (central charge c = 1 for the 
corresponding IR CFT.), and the results support the scaling law (fl~3l) and the entropy bound 

(ED. 

As mentioned above, the straightforward calculation of block entanglement entropy is 
formidable because of the exponential growth of the block's physical dimension. Instead, 



we follow the method of quantum state RG transformation proposed in [3CJ so that the 
physical dimension is capped by x 2 , be., the correlation length of the MPS, and then we can 
obtain the block entanglement entropy of size L by performing L steps of quantum state RG 
transformation. However, we will see that the quantum state RG transformation method 
fails to get the expected area law for the block entanglement entropy in 2D case due the 
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exponential growth of bond dimension when merging sites. 
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FIG. 9: The GE of block size L for spin 1 XXZ chain with the fixed J z = 2.59 and the quantum 
critical point at D ~ 2.3. We adopt the MPSs of x = 2, 4 for the ground state. (Left) Off the 
critical point, the block GE saturates around L = 4. (Right) Very near the critical point, the block 
GE saturates around L ~ x 2 , which can be seen more clearly from the numerical data. 



On the other hand, we like to see if the other entanglement measures will also obey 
the scaling law CLU]) and the entropy bound (fT4"l) . We then evaluate the block GE for the 
spin 1 XXZ chain. Again, the direct computation of block GE is formidable because of 
the exponential growth of physical dimension for the trial separable states. Instead, one 
can adopt the method of quantum state RG transformation. In this way, one can still 
use one-site trial separable states but need to optimize the two-site unitary transformation 
for the quantum state RG transformation. This optimization only involves a vector of x 
components, and it is numerically tractable if x is n °t very large [4^]. Our result for the 
block GE of the spin 1 XXZ chain with x = 2, 4 is plotted in Fig. 9. Again, we find that 
the results support the scaling law ( 1T31) and the entropy bound (fT4"l) . 

For comparison, in Fig. 10 we also plot the bipartite entanglement as a function of x, we 
see that it does saturate off critical point, but the scaling law is £ x ~ ^- 2 - 286 instead of ffT3l) . 
It is not clear if the deviation is due to the truncation error of the approximate MPS or not. 
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FIG. 10: The scaling behavior of bipartite entanglement for spin 1 XXZ chain. The critical point 
is at D ~ 2.3. 

B. 2D scaling via TPSs ? 

We now aim to see if we can find the scaling law of the entanglement measure from the 
TPSs for 2D quantum spin systems. The direct computation cannot work as in the ID case, 
and we then wonder if the quantum state RG transformation works for our purpose or not. 

As in ID case, to keep the physical dimension bounded after merging four neighboring 
sites into a new block by TRG, we then need to perform some appropriate SVD for the 
resulting tensor to obtain the non-local unitary transformation. Its physical dimension is 
bounded by x 4 - However, unlike the ID case, the bond dimension will scale up exponentially 
after each merging, it implies that the bound on physical dimension will scale up accordingly. 
In fact, this is the same reason to introduce TRG method to calculate the double tensor 
trace by truncating the bond dimension, except now we do not contract out the physical 
indices when the performing quantum state RG transformations. 

To be more specific, let A s rdlu denote the tensor of a site for TPS with s the physical 
index, and r,d,l,u the bond ones. After merging by usual TRG method, the resulting block 
tensor is denoted by A^'**'* 3 . Note that (s\, S2, S3, S4) is the combined physical index 
whose dimension is d*, and the bond dimension is fixed as in the usual TRG method so that 
the bond indices r,d,l,u remain intact. Then, we perform the SVD so that A^'**' 83 ' 34 ^ = 
^mm(d , x ) jjKsi,s 2 ,s3,s4,) \ s y& ^ to bound the physical dimension by x 4 , or more precisely by 
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mm(dg, x 4 )- After that, the non-local unitary operator JJ^ sl ' S2 ' S3 ' Si ' ) is used to RG transform 
the quantum state. Finally, we need to perform TRG again to keep x fixed when evaluating 
the double tensor trace for the entanglement measure. 





FIG. 11: Spectra of singular values for initial bond dimension x = 4- The index a is to number the 
singular values Aq,'s. (Left) The spectrum of singular values for bond dimension after performing 
the TRG to merge the four neighboring sites. (Right) The spectrum of singular values for physical 
dimension after performing the SVD to find the non-local unitary transformation for quantum 
state RG. 

Follow this prescription, we compute the block von Neumann entanglement entropy of 
the 2D transverse Ising model up to x = 4. However, we find that the resulting block 
entanglement decreases as the block size L increases. We have also found that the resulting 
value of entanglement is not so sensitive to the value of physical dimension. This failure 
may mean that we truncate the bond dimension too much in order to keep x fixed when 
performing quantum state RG transformation, so that the entanglement is lost even when 
we increase the block size. Of course, we could try to increase the x value to keep more 
entanglement intact, but the computational cost for doing this is beyond what we can afford 
at this moment. Instead, in Fig. 11 we plot the spectra of singular values after one step 
of quantum state RG transformation, which can indicate the loss of quantum entanglement 
due to the truncation of the bond and physical dimensions. It shows that we want to keep 
the singular values greater than 10~ 3 , we should fix x — 10 which, however, is beyond our 
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computation power. 



V. CONCLUSION 

We have studied the global entanglement measures and bipartite entanglement entropy 
for various ID and 2D quantum spin systems. This is achieved by the iTEBD and TRG 
algorithms. Like the 1-tangle or concurrence, we find that both the global entanglement 
measure and the bipartite entanglement entropy are capable to characterize the quantum 
critical points as good as the order parameters. Moreover, they sometimes can capture the 
non-entangled states which are missed by the conventional order parameters. Especially, 
the bipartite entanglement entropy is far easier to compute in the MPSs/TPSs formalism 
than the most of the other entanglement measures. Interestingly, the bipartite entanglement 
entropy is always larger than the global entanglement measure for the cases studied here. It 
is hoped that one can understand if there is a universal class for the quantum entanglement 
measures, and the relation among the members of such a class. 

Besides, we also study the scaling behaviors of entanglement measures near the quantum 
critical point, which again shows some universal feature at least for the ID spin systems. 
However, the method of the quantum state RG transformation fails to capture either the area 
law or the scaling law of the entanglement measures for 2D systems. It deserves more works 
to elucidate the subtlety of truncating the bond and physical dimensions, and we hope to 
find a viable numerical algorithm in carrying out the scaling behaviors for 2D entanglement 
measures. 
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